Quantum discord of two-qubit X-states 
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Quantum discord provides a measure for quantifying quantum correlations beyond entanglement 
and is very hard to compute even for two-qubit states because of the minimization over all possible 
measurements. Recently a simple algorithm to evaluate the quantum discord for two-qubit X-states 
is proposed by Ali, Rau and Alber [Phys. Rev. A 81, 042105 (2010)] with minimization taken 
over only a few cases. Here we shall at first identify a class of X-states, whose quantum discord 
can be evaluated analytically without any minimization, for which their algorithm is valid, and also 
identify a family of X-states for which their algorithm fails. And then we demonstrate that this 
special family of X-states provides furthermore an explicit example for the inequivalence between 
the minimization over positive operator- valued measures and that over von Neumann measurements. 
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It is believed that entanglement is an essential resource 
in quantum computation and communication. However, 
there are also quantum tasks that display the quantum 
advantage without entanglement, for example, the deter- 
ministic quantum computation with one qubit [l| . There- 
fore, there is a need to consider quantum correlations be- 
yond entanglement 0]. The quantum discord has been 
introduced by OUivier and Zurek Q and independently 
by Henderson and Vedral Q to quantify quantum cor- 
relations. Recently the quantum discord has attracted 
much interest in quantum information theory [H-Q such 
as its relation with the complete positivity [lO| and local 
broadcasting of the state [11[. Furthermore, both Marko- 
vian and non-Markovian dynamics of quantum discord 
have been analyzed not only in theory but also in ex- 
periments |12Hlq . Whether the quantum discord can be 
more robust against decohercnce [3 than entangle- 
ment or not [l7j is still an open problem [l7| . 

The quantum discord is always nonnegative Q . States 
with vanishing quantum discord are relatively well un- 
derstood and necessary and sufficient conditions are ob- 
tained to detect nonzero quantum discord [l^ as well as 
nonlinear witnesses have been proposed both for a given 
state and an unknown state Unfortunately al- 
most all quantum states have nonzero quantum discords 
p^ . which arc notoriously difficult to compute because 
of the minimization over all possible positive operator- 
valued measures (POVMs) or von Neumann measure- 
ments. In addition to a few analytical results including 
the Bell-diagonal states [2l| , rank-2 states [l^j , and gaus- 
sian states [2^ , a thorough numerical calculation [2J] has 
also been carried out in the case of von Neumann mea- 



surements for two-qubit states. 

For an important family of two-qubit states, the so 
called X-states [1^, an algorithm has been proposed to 
calculate their quantum discord with minimization taken 
over only a few simple cases 12611 . which is unfortunately 
impeded by a counter example |27| . In this paper we shall 
at first identify a vast class of X-states, whose quantum 
discord can be evaluated analytically without any min- 
imization at all, for which their algorithm is valid, and 
also identify a family of X-states Xm, the so-called max- 
imally discordant mixed states [l^l, for which the above 
mentioned algorithm fails. And then for this family of X- 
states Xm we construct a POVM showing that the quan- 
tum discord obtained by minimization over all POVMs is 
strictly smaller than that over all possible von Neumann 
measurements . 

For a given quantum state p of a composite system 
AB the total amount of correlations, including classical 
and quantum correlations, is quantified by the quantum 
mutual information T{p) = S^qa) + S{gB) — S{g) where 
S{g) = — Tr(gilog2 g) denotes the von Neumann entropy 
and gA, gs arc reduced density matrices for subsystem 
A, B respectively. An alternative version of the mutual 
information can be defined as 



Me) 



Sign) - min ^PkS{gB\k) 



(1) 



where the minimum is taken over all possible POVMs 
{E^} on subsystem A with pk = Tr(i?^p) and gB\k = 
TYA{E^g)/pk- Since jA{g) quantifies the classical corre- 
lation, the difference 



DAig) = i{p)-Ja{q) 



(2) 
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defines the quantum discord that quantifies the quan- 
tum correlation. Also the minimum in Eq.ffl can be 
taken over all von Neumann measurements |3| and we 
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denote the corresponding classical correlation as Ja{q) 
and quantum discord as Da{q), respectively. Obviously 
Da{q) < Da{q) and it becomes an equality for some 
states such as Bell-diagonal states and a family of filtered 
X-states [111. 

The two-qubit X-state usually arises as the two- 
particle reduced density matrix in many physical systems 
possessing z-axis symmetry. In the computational basis 
{|00), |01), |10), 111)} its density matrix 



X = 



goo gos \ 

Qii Qi2 
922 



(3) 



has seven real parameters. Via local unitary transforma- 
tions, which preserve the quantum discord, elements qq^ 
and Q12 can be brought into real numbers. Thus there 
are in fact only five real parameters, which can be con- 
veniently taken as 



X = Qm + Qii - Q22 - Q33 = Tr(cr2 X), 

V = Qm- Qii + Q22 - 933 = Tr((Tf X), 

t = Qm- Qii - Q22 + Q33 = Tr((T^CTf X), 

s = 2(f?i2 + te) =Tr(a^afA'), 

u = 2(^?l2-^?03) =Tr(a^t7^A'), 



(4) 



where cr^'-^z three standard Pauli matrices. The pos- 
itivity requires (1±<)^ > (x ± y)^ + {sT u)"^ with all five 
parameters taking values in the interval [—1, 1]. Without 
loss of generality we shall assume \s\ > \u\ in what fol- 
lows because we can always change the sign of go3 by a 
local unitary transformation. 

The quantum discords of some X-states have been nu- 
merically calculated [13, [IB 13 ■ Most recently an algo- 
rithm j26| has been proposed to calculate quantum dis- 
cord for all two-qubit X-states in which the minimization 
is taken over only a few simple cases instead of actually 
minimizing over all possible measurements. In the case 
of real X-state with \s\ > \u\ the algorithm reads: for 
Da{X) 



the optimal observable is either or cr^. 



(5) 



Recently a counter example is found in [27| for this algo- 
rithm. The following theorem identifies a region of pa- 
rameters of X-state, whose quantum discord can be eval- 
uated analytically without any minimization, for which 
the above mentioned algorithm is valid. 



Theorem The optimal measurement for the quantum 
discord Da{X) and Da{X) of a real X-state X with \s\ > 
\u\ is i) if 



(l£'i2| + \Qa3\)^ < {goo - fiiiOCte - £'22) 



(6) 



and ii) if |^ goo £'33 - ^/QlTQ22\ < \gi2\ + Ifoal- 

Proof Considering a general POVM {/ife(l + 
^k^^)}^=i with K < A [s^l made on the first qubit, 
where J^k = 1^ "fc = 1: ^^d J^k f^kfik = 0, we obtain 
the outcome k with probability pk = ^k{^ + xukz)- The 
second qubit is in the conditioned state 



'^Blk — 



1 + xrikz + snkxcr^ + uukya^ + {v + nkzt)cr^ 
2(1 -I- xukz) 



•(7) 



By denoting h{w) = loga ^ - ^ loga ^ and 
Afc = (1 - nlJs'^ + {y + rikztY, we obtain 



Pkh. 

k k 



flkV^^ 



s. 



Pk 



X, 



(8) 



where the inequality is due to the function h{'w) being 
a decreasing function for w > and the equality can 
be attained by taking Uky = 0,Vfc. Note that in the 
case of von Neumann measurements, since any observable 
(T^ not in the x-z plane the observable a^, with n' = 
(•^1 — rt^, 0, riz) will yield a smaller value, the optimal 
von Neumann measurement must lie in the x-z plane. 

Note that in the case of von Neumann measurements, 
any observable not in the x-z plane the observable 
^n' with n' = (-^1 — n^, 0, riz) will yield a smaller value. 
Thus the optimal von Neumann measurement must lie in 
the x-z plane. 

By denoting Afe± = /-tfc(l + xukz ± VAfc)/2 and its 
marginal Afe = Xk+ + Xk- = Mfc(l + xukz), we obtain 



Sx =^ Afc log2 Afc - ^ Xk± log2 Afc± . 



(9) 



k± 



In the following we will find out the minimum of Sx for 
fixed {^1, ...,^k}- Since X^fcMfc = 1 and Mfe ^ we can 
suppose without loss of generality that fiK > 0. Be- 
cause of condition J2k ^^k^^'kz = we can regard ukz as a 
function of if — 1 independent variables nkz '■= nk with 
A: = 1,2,..., A'— 1. Therefore 5*;^ is a multivariable func- 
tion of Tifc with k ~ I, K — 1 and its Hessian matrix. 



whose elements are VFi< 



driidn-i 



reads 
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K 



E 



k=l 



fc=l,± 



E 

fc=i 



AfcAfc+Afc- 



fe± 



y2 _ ,2) 



In 



(10) 



where A 



k± 



k± \i] 



and Al = \k-X'k+ - Afe+Aj,_. 



driidnj ' 



Bx(m,e) < ^ 
J^^^ Bz(m,e) < 



1 

10 



3 
10 



FIG. 1: (Color online) The optimal observable for the state 
parameterized by m and e is i) (Jx if Bx (jn, e) < (top 
red region) and ii) if Bz{m,e) < (bottom blue region). 
The black curve corresponds the state Xm- 



If ^ + which is equivalent to Eq.®, the matrix 
W is always negative semidefinite, thus the conditional 
entropy Sx is a concave function of riiz- Therefore the 
minimum of Sx is attained on the boundary, i.e., tiiz = 
— 1 or riiz = l,Vi. Thus for every given {^i, the 
optimal measurement is a^, which proves the first case 
of the Theorem. For the proof of the second case, we 
refer to Hg. Q.E.D. 

As the first example we consider the Bcll-diagonal 
states for which we have x = y = 0. If |i| > \s\ then 
we have case i) so that the optimal observable is cr^ and 
if |<| < \s\ then we have case ii) so that the optimal ob- 
servable is ax (recalling that we have assumed \s\ > |m|), 
both reproduce the result in [IJ] . As the second example 
we consider a 2-parametcr family of X-states 



/ e/2 

(l-e)m 

(l-e)(l 

\ e/2 



(11) 



with < e, 2m < 1. In this case we have x = —y = 
(1 — e)(2m — 1), s = —u = e, and t = 2e — 1. According 
to two cases in the Theorem the optimal observable is 



0,229 



I" 0,228 



0,227 




0.102 



FIG. 2: (Color online) The optimal observable axSinOopt + 
ffz cos opt for the quantum discord Da{X3). 



ax or CTz in the case of Bx{m,e) < or Bz{m,e) < 0, 
respectively, where 



Bx{m,e) = ^m(l - m) 



Bz{m,e) 



1 - e 



1-e' 
2m(l — to). 



(12) 
(13) 



If e > 1/3 we have Bx{m, e) < and if e < 1/3 we have 
illustrated in FigUthose two regions (top red and bottom 
blue) of parameters for which the optimal observable is 
either ax or az and in these regions the algorithm ([SJ is 
valid. 

However there also exists a region (white) of parame- 
ters about which our theorem docs not say anything. To 
find out the quantum discord £)a('%3) in this region we 
have to do numerical calculations. Taking into considera- 
tion the eigenvalues A (A'a) = {e, (1 — e)TO, (1 — e)(l— to), 0} 
and XlTrsX^) = {(1 ± x)/2} we obtain the value 

De = Sxs (cos 0) + h{x) - h{t) - (1 - e)h{2m - 1) (14) 

after the measurement of the observable along an arbi- 
trary direction n ~ (sin0, 0, cos 0) in the x-z plane. Here 
Sx^ (cos 0) is defined in Eq.® with riz = cos 0. The mini- 
mization of De over all possible angles gives the quantum 
discord, i.e., Da[Xz) = ming Dg Dg^^^. We note that 

Do = e and SxM = /i(VF+^)- 
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Doapt (dashed red line), and D^^"^^ [Xm) (dotted red line). 



FIG. 3: (Color online) As the functions of e, the optimal 
angle 6 opt for Z),4(A'm), the difference A = e — DA{Xm), and 
the difference A = e — D^^'^^ [Xm) are plotted in solid blue, 
dashed red and dotted red curves, respectively. 



A detailed numerical search for DAi^^) in the ranges 
m e [0.1,0.102] and e e [0.227,0.229] has been car- 
ried out with results shown in Fig. 2. The parameters 
for which the optimal observable is ax or ct^ are high- 
lighted in red (top) or blue (bottom) respectively. How- 
ever these two regions are separated by an intermediate 
region for which the optimal observable is neither Gx nor 
CTj. As a result, though valid for most of the parame- 
ters, the algorithm ([S]) fails for a region with finite mea- 
sure. Moreover our numerical search shows that there 
are about 0.05% of 5 x 10^ randomly chosen X-states 
satisfying s ^ u that violate the algorithm ([S]). In our 
search the difference between the true value of Da{X) 
and min{Z?Oj D^} can reach as high as 0.0029 in the case 
of x = -0.8812,'y = 0.9407, s = 0.2898, t = -0.9383. 

Especially we have considered in more details a sub- 
family of X-state Xm of X^^ with e determined by the 
condition Do = oi more explicitly 



h{^Jx'^+e^) + h{x) - h{t) - (1 - e)h{2m - 1) = e (15) 

recalling that a; = (1 — e){2m — 1) and t ~ 2e ~ 1. This 
special family of X-states Xm is exactly a familyof so- 
called maximal discordant states investigated in [24| . As 
a function of m the solution e{m) of the above equation 
is plotted as a black curve in Fig.l. 

In FiglH the optimal angle 9opt (solid blue line) and 
the difference A = Z^o ~ ^Sopt (dashed red line) are plot- 
ted as functions of e. As e{m) increases, A increases 
at first from to its maximum about 1.07 x 10~^ at 
e ~ 0.115699 and then decreases to zero. Except at the 
end points (to, e) = (0, 0), (1/2, 1/3) for which any ob- 
servable is optimal, we have A > which means that 



and arc not optimal. Most interestingly the optimal 
angle 9 opt takes values in [0,7r/4] continuously. This fact 
strengthens the theorem proposed in Ref. j27| . i.e., it is 
impossible to find a universal finite set of optimal von 
Neumann measurements even for the real X-states. 

In Ref. [2^ . by maximizing the quantum discord for 
given classical correlations over all states A3, whose quan- 
tum discord is taken to be min{Do,^f }, the state Xm 
turns out to be a family of maximally discordant mixed 
state that lies on the boundary of the J"-D diagram 
of the classical-correlation vs quantum-discord. Since 
-D^(A'm) < uim{Do, D^} except at the endpoints, the 
J^-D diagram for Xm is shifted right-downward as shown 
in Fig. m Thus whether the state Xm is still on the 
boundary or not needs further numerical calculation to 
substantiate. 

Finally let us consider the quantum discord DA{Xm) 
obtained by minimization over all possible POVMs. For 



each given state X„i there is an optimal angle 



^opt 



for 



DAi'^m) obtained by von Neumann measurements and 
we perform a 3-outcome POVM {Mfe(l+"'fc'?'^)}fe=i made 
on the first qubit on the subsystem A, where 



cos 6, 



1 + cos 9opt 
1 



, ni= {0,0,-1}, 
, n2 = {sin 9 opt, 0, cos 9 opt}, (16) 



2(1 + cos 6'opt) 

^3 = TTTT—-— 5 — T' na = {-sin 6'oj,t,0, cos 6'opt}. 
2(1 + cos 6'opt) 

We denote by D^^'^^ [Xm) the corresponding suboptimal 



value for the quantum discord and obviously DA{Xm) < 
^7^"'(A'™). The difference K ^ - Dj^'"' {X^) is 
shown in Fig. [3] as a dotted red line and we have A < A 
(except at end points), which means that the S-D dia- 
gram for Xm must be shifted further right-downward as 
shown in Fig. 0] by the dotted red line. The boundaries 
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of the J-D diagrams for POVMs and von Neumann mea- 
surements would be different if were still the maxi- 
mally discordant states. 

To summarize, we have presented some positive results 
as well as negative results on the quantum discord of X- 
state. We have identified a vast class of X-states whose 
quantum discords can be evaluated analytically and also 
a family of maximally discordant mixed states Xm that 
invalidate the algorithm [2^ . If the state X^ were still on 



the boundary of the J-D diagram [2J| then this part of 



boundary would be shifted right-downward and even fur- 
ther for POVMs. Thus the state Xm provides an explicit 
example for the inequivalence between the minimization 
over POVMs and that over von Neumann measurements 
for X-states. Recently, more examples [3l| have been 
given. 
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